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SUMMARY 


The effects of diffusion damping of spin-echo amplitudes 
were first treated by Hahn’, An improvement of Hahn's 
averaging procedure was suggested by Das and Saha“. Their 
theory, however, as is shown here, contains an error. This 
is corrected, and the revised theory is shown to give for a 
Special case the same result as obtained by different methods 
by Carr and Purcell? and by Torrey’. Furthermore the 
revised theory is extended to account for any combination of 
pulses, and it is shown that the results are conveniently 


summarized in a simple "spin-echo diagram". 
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GENERAL THEORY ON SPIN- ECHOES FOR ANY 


COMBINATION OF ANY NUMBER OF PULSES, 


INTRODUCTION OF A SIMPLE "SPIN-ECHO DIAGRAM" 


Erling Jensen 


Technical University of Denmark 


Introduction 


The damping effects of diffusion on spin-echo amplitudes were 
treated in Hahn's original work’, However, an improvement in Hahn's 
averaging procedure was suggested by Das and Saha” who pointed out 
that the damping due to diffusion after a single pulse might be affected 
by successive pulses, The special case ofa 90° pulse followed by a 
sequence of 180° pulses was treated by Carr and Purcell? who used a 
simple and convincing random-walk theory. Application of the theory 
of Das and Saha on the Carr-Purcell pulse method, however, yields 
a result which is not in agreement with that obtained by Carr and 
Purcell. As will be pointed out here the source of this discrepancy 
is that in the theory of Das and Saha two sets of statistical variables 
are treated as if they were mutually independent when, in fact, they 
are directly correlated, 

Here will be given a theory, which takes account of this correla- 
tion and may be characterized as a combination of the theories of 
Das and Saha and of Carr and Purcell. It gives as a special case the 
resultof Carrand Purcell but furthermore is used to derive a general 
solution for all the spin-echoes due to any combination of pulses. 

The results are most conveniently summarized in a simple "'spin+recho 


diagram" that will be introduced later on. 


— 


Summary of previous theories and critical discussion of same 


For later use it will be convenient to notice the following relations 


and to summarize some previous results. Let 


n 
x=) x; (1) 
is1 


where the X; are mutually independent statistical variables with even 


density functions p(X;). Then 


co 


n n 
E(cos X) -{ cosX p(X) dX = qT E(cosX.) = TT { cosX. p(X.) dx; (2) 
as all other terms contain at least one independent factor sinX,, and 
thus integrate to zero, If furthermore the X; are normal random 
variables 


n 


E(cos X) = exp |- 5 E(x’) = exp |- ry) E(x,”) ; (3) 


Only taking account of self-diffusion with the self-diffusion coeffi- 
dH 
cient D and assuming a constant field gradient G, (G = ar) 


= YGl 4 
Y tt, (4) 


E(n, , = - t,) = 2k(t - t,) (5) 


et is the diffusion length in the 


time interval t - tye and "eet is the corresponding random increase 
1 


where y is the gyromagnetic ratio, 1 


in angular velocity. 
Then the random phase accumulation at time t after a single pulse 
supplied at time 0 is (see ref. 3) 
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Since the n, are mutually independent statistical variables with 


mean value zero, and by using (5) 
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Fig. 1. Representation for 


For simplicity the array (6) for , ™ will merely be symbolized by 
the triangle fig. 1 containing the n; and characterized by the indicated 


time axis. 
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Fig. 2. Applied pulses at the radio-frequency. 


The pulse array under consideration is shown in fig. 2, and in 
table I (obtained by slightly modifying table I, ref. 2, pag.751, which is 
derived directly from-the Bloch equations) are given the normalized 
components U, V, and W of the magnetization vector M(Aw) when the 
values are averaged with respect to Aw and suitable diffusion variables. 
To express the random part of the phases hes ty Das and Saha make 
use of both 6's and n's (the first form at the bottom of table I), whereas 
it is preferred in this paper only to use 's (the second form). 

Assuming Aw to have an even density function the U-component 
integrates to zero as it is seen always to be an odd function of Aw, which 
is of course independent of any diffusion variable. The signal therefore 
is simply V(t). Thus the free induction decay following the first pulse, 
Oo<t< Ty» is found from table I replacing T, with t, and according to 
the above characteristics of Aw using (2), and thereupon (3) and (7), 
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Assuming here for simplicity Aw to be normally distributed with 
variance ) » (3) gives directly 
3 


t 1,3 
~sin£ ,exp(- (9) 
3 


For the primary echo Hahn now simply uses the diffusion factor 
for the decay following a single pulse, inserting t = 27, 


8 3 
o) = exp(- ) (10) 
Das and Saha, however, assuming T and be 


mutually independent, even normally distzibuted statistical 
obtain for the primary echo, and for the (by the second pulse modified) 


free induction decay from the first pulse, the diffusion factor 
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Component- sums of >, i after Das and Saha. 
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Considering the array (6) for o T 


(see fig. 3) and from (3) it then immediately follows that Das and Saha 
would arrive at exactly the same result as Hahn if a oe *ey,, 7, and 
Nr. 9 Were mutually independent, that is if the condition for their 
calculations did hold. Since in the array for >, - the same n is repeated 
in vertical direction while y's in a horizontal line are mutually independent 


is of course independent of both , ‘~ and of however, 
4’ 


o, and Nr. 9 are directly correlated. (In general, any two n or 
49 4° 


are always correlated when they belong to time intervals overlapping 
each other; otherwise always uncorrelated. In other words the triangle 
for >, o May be divided into mutually independent component-sums by 


vertical, and only by vertical, lines). 


Corrected calculations 


In the following will be given the revised theory taking account of 
the possibility of correlation between component-sums of >, ” From 
table I the signal V(E) for T, <t <T, >27,, replacing T, witht, is 
found to be: 


exp(- 


“cose, cosA_ cosA, 
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The bracket from the first term gives 
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Thus V(E) is: 
252 t 
~sin— ,cos exp (- [eos(aut + 


exp (- DIE - 2T,) + (o, 26, | 
2 4° 
t- 


or, by using (3) and regarding fig. 4 indicating how the resulting diffusion 
sum — (enclosed in the heavy lines) may be divided into mutually indepe«- 


dent component-sums 


V(E) = 
~sin£ ,cos exp(- —;)exp(- (12a) 
2 2T3 
(t-27,)° 1 
+eink , sin’ exp (42b) 
t-T (t-T,) 
-W(T, )sint , exp(- — exp (12c) 
3 


E 
(12a) is the fraction cos 5 of the free induction from the first pulse 


which continues with unaltered diffusion term. 


(12b) is the primary echo formed by the second pulse from the rest of 
the decay from the first pulse. The echo has reduced diffusion 
damping, and in the special case of echo maximum t=2T, the factor 


is exp(- Sur, *) in agreement with Carr and Purcell and with Torrey’. 
(12c) is the free-induction decay following the second pulse. 


As the normal spin-echo conditions are of course assumed to be 
satisfied, inserting Eo = 360° is equivalent to omitting the second pulse, 
and as one would expect, the total free induction decay from the first 
pulse is seen to continue unchanged, while the echo with reduced diffu- 


sion damping does not occur. Das and Saha, however, obtain for both 
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Fig. 4. Illustration of — for the three components of the 


signal following the second pulse. 


the echo and the modified decay the same reduction of diffusion damping 


(see (141)) which does not seem reasonable. 


Performing quite similar calculations V(H), T, < t< T3 > 2T>, is 


found to be: 
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where 


T,°T 


Again by using (3) and fig. 5 for j = 3 (an extension of fig. 4) 


V(H) = 
E 2 
-sint exp(- exp(- exp(- Sut?) (13a) 
3 
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2 
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x exp + (t-T,47, } (43h) 
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(t-T,-7 
+4eint sin’ ,exp(- jexp |- 
2 2T 
3 
k 3 2 3 , 
x exp| - 3} +3(t-7,) (72-74) - (t-72-7,) }] (413i) 


Here (13a, b, e) are modifications of (12a, b, c) with unaltered diffusion 
terms, (13c,d,f) are echoes of (12a, b, c) caused by the third pulse, (13g) 
is the free-induction decay following the third pulse, (13h) is a virtual 
and (13i) a real stimulated echo. Table II gives for the four secondary 
echoes the diffusion terms at echo maximum (compare table III, ref. 2, 
pag. 754). 


Table II. Corrected values of diffusion terms for the secondary echoes. 


Echo Value of diffusion term 
Term maximum | at echo maximum 
Stimulated echo ex [- K(37%7 ~ 73)| (413i) 


First 

secondary echo} exp + (7, - (13d) 
Second 

secondary echo | 27,-T, exp | - (13) 
Third 

secondary echo | 27, exp | (13c) 


It is now obvious how to continue with four pulses or more by 
extending table I. On account of the striking systematic nature of the 
results, however, it will be convenient to proceed directly to the 
following recursion formula, obtained simply by analysing the "transport" 


of terms in the extended table I. For * t< 
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and where B. ,=B. ,(T.), starting with: 


j 
T, 
Bi By. 47 By 4 (74) = ,exp(- T,) = , exp(- T,) 
The number x(j) of signal components for <t<T,> 27 4,2, 3..«) 


is seen to start with: x(i)= 1; x(2) = 3; etc. Likewise, 
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Fig. 5. "Spin-echo diagram" (for four pulses or less). 
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(14) leads to the simple ''spin-echo-diagram"' fig. 5 from which any 
information about any signal component for a wanted number of pulses may 
be obtained according to the following directions. Starting as indicated for 
a wanted pulse number j and choosing, when crossing each of the vertical 
lines corresponding to pulses, one of the three possibilities, (the horizontal 
or one of the two inclined lines), the ''maximum times"! ts . for the number 
x(j) of signal terms Ti04 <t< " are simply found as the points where the 
system of broken lines crosses the time axis (which is the vertical line 
corresponding to the first pulse). Thus the total number of component 


terms due ton pulses is: 


n n 
j=1 j=1 
However, x(j) is made up of real echoes, virtual echoes, and of pulse 
decays. If T.> 7 gs as assumed for simplicity in the previous calcula- 
tions and for fig. 5, the total number of real echoes is: 
n n j n 
k=2 1,.j-1 
N(a) = y@) » 4-4) = $(3"-2n-4) (17) 
j=2 j=2 k=2 j=2 


which result has been derived before (see ref. 5). 

The trigonometrical coefficients of the different signal terms are also 
obtainable from the diagram, observing the following rules: Immediately 
after 'start'' in the diagram, and also whenever moving along an inclined 


line, then upon crossing the j'th vertical pulse line, 
a) if continuing in the same direction: multiply by cos” bs 
b) if choosing the other inclined line: multiply by @sin” = 
c) if choosing the horizontal line: multiply by ~sing ; x 


When coming from a horizontal line and crossing the j'th vertical 


pulse line, 
a) if continuing in the same direction: multiply by coet, 


b) if choosing one of the inclined lines: multiply by seing, : 


However, if after crossing the j'th vertical pulse line the path 
continues in a horizontal direction, it is necessary to insert the factor 
~sinE, should be inserted. The K's are the only factors of the expressions 


Especially by intersecting the time axis the factor ~K 


for the various signal components which are combined trigonometrical and 
(longitudinal) relaxational terms (see (15)). 

Pure exponential transversal relaxation damping is present during 
the total time corresponding to inclined lines. Pure exponential longitu- 
dinal relaxation damping is present during the total time corresponding to 
horizontal lines, (excepting a horizontal ''K-line'' touching the time axis, 
for which as mentioned the longitudinal damping is not simply exponential). 


Only the diffusion damping remains to be considered. The factor for 


this is seen to be exp E | » where j res(t) = o j 

is the resulting random diffusion sum enclosed between the horizontal line 
which corresponds to the wanted time t, and the broken line which corresponds 
to the considered i'th number of signal components following the j'th pulse. 
The geometrical figure containing a is divided by vertical lines into 
rectangles, triangles, and trapezoids containing mutually independent 
normally distributed diffusion sums. As the variance of the sum within 

a trapezoid is the difference between the variances of the sums within two 
triangles, the diffusion factor is easily obtained by using (3) together with 

(5) and (7). 

As a correction of previous results obtained by Das and Roy (see ref. 6, 
table I, pag. 526) the signal components following the fourth pulse, 
T3<t<T,> 2T 3, will now be taken directly from the spin-echo diagram 
fig. 5 (The more or less modified decays from previous pulses and echoes 


are omitted as they are of very little interest), 
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So far it has been assumed that T. > aT. 4s but from the very way in 
which the calculations are performed it follows that t , and $5 tte (t) 
will be given by the same expressions, whichever pulse train is used. 
So also a spin-echo diagram like fig. 5 remains valid for any combinations 
of pulses, and the expressions for the various terms can still be derived 
from the diagram, following the above given directions. The only difference 
resulting from supplying new pulses before all the echoes from the previous 
one have formed is that the number N(n) of real echoes is reduced (as some 
of the otherwise real echoes now become virtual), and that the echoes formed 
become further (triginometrically) modified. 

For an echo which is the result of repeated echo formations it will be 
seen that the diffusion damping at echo maximum is the accumulation of 
the diffusion dampings for each echo formation. This follows simply from 
the fact that — is enclosed between the broken line for the considered 
echo and the horizontal time-line at echo maximum, and echoes occur 
when those two lines are crossing each other. 

For the Carr-Purcell pulse train ,=90°, gud, 3,204, 8) 


the trigonometrical part of the coefficient is: 
,(-1) sin = = (-4)p. 


Theoretically it is possible to adjust to this value if the pulse train is 
assumed to be stationary, and if the radiofrequency field is sufficiently 
homogeneous. If, however, inhomogenity is present the ideal rotation 
angle cannot be obtained for all of the spins, and a damping (besides 


possible adjustment damping) will occur, the factor being: 
j 


1-E(AE , for j>>1 and AE. = AE. 


By means of a spin-echo diagram, however, it is easy to see that with 
the ratios of pulse intervals, 1:2:2:2:2 etc. used by Carr and Purcell 

a lot of disturbing echoes willoccur exactly at the maximum times for 
the considered echo array. Among the disturbing echoes one group is 
quite dominating, namely stimulated echoes (and repetitions of same) 
caused by the pulse train together with the first pulse and the considered 
echo array. They will all have the size ~ FE (at), and will all have 


20 


the same sign which is opposite to but shifting in accordance to that of an 


echo of the considered array. Since the number of disturbing echoes is 
$(-4) (j-2) Si for j>>1 they will give a damping factor: 
Thus the fault due to disturbing echoes largely exceeds the fault on the 
considered array itself. (Since T, 2 T, this effect is normally increased 


by taking account of relaxation). It may be of some interest to note that 


for the same pulse train modified to have the pulse interval ratios 


2:6:7:7:7 etc. all disturbing echoes from a spin-echo diagram are seen 


to be separated from the considered array. Furthermore their size may 


be used to some extent to indicate whether the pulse angles have the ideal 


values or not. 
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